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Abstract 

We generalize Schwinger's original mass formula to the case of an additional 
isosinglet mixing with the nonet mesons. We then make further generalization 
to either (i) an arbitrary number of additional isosinglets mixing with nonet 
mesons, or (ii) an arbitrary number of mesons mixing with an additional isos- 
inglet. In the former case, we present an explicit relation, while in the latter 
case, we show by numerical example that the new mass formula is only weakly 
affected by the inclusion of additional mesons, and hence holds with good ac- 
curacy for each of the 3x3 mass sub-matrices. 
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1 Introduction 



Schwinger's original mass formula [||] has attained celebrity due to its generality in 
a universe with only nine mesons, under the assumption that up, down and strange 
quarks may be treated the same. 

If there were only up and down quarks which interact identically via the strong 
force, the uu and (id states of the hadron spectrum could be resolved into degenerate 
isovector and isosinglet statesQ Although this requires more than isospin invariance, 
we know it is consistent (for mesons) with the details of QCD dynamics as there are 
no spin-isospin correlations here to produce color magnetic spin interaction induced 
splittings, at least for the lowest Fock components. 

The isovector masses can be regarded as physical since they undergo no further 
mixing. Therefore, in the Schwinger formula, the primitive (bare) isosinglet masses 
can be equated with the (physical) isovector masses. Isosinglet states are invariant un- 
der rotation of up into down quarks and have the flavour structure nn = (uu+dd) / \[2. 
Enter the strange quark: Primitive ss states, the masses of which are not a priori 
known, can mix with the primitive nn states to create physical states. Schwinger's 
formula pertains to the case where there is one nn and one ss primitive state, and 
where the strange quark behaves the same as the up and down quarks when the two 
primitive states couple to each other (SU(3) symmetry). 

The formula has the benefit, except for the primitive ss mass, of referring only to 
experimentally measurable meson masses, and not intractable couplings. 

In this Letter we extend Schwinger's formula to include additional isosinglets and 
relax the SU(3) assumption of Schwinger. This enables a more model independent 
analysis of the presence of isosinglets. In the case where the isosinglet is built mainly 
from gluons (a glueball), this extension is of central significance as some experimen- 
tally discovered isosinglet zero angular momentum states are very widely accepted 
to exhibit mixing with a glueball |3|, [|]. 

Schwinger's formula 



3M(ss)-2M(ss)-M(nn) 3M(nn) -2M(ss) - M{nn) +2 M{ss)-M{nn 
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o. 

(1) 

relates the masses of the physical (M(ss), M(nn)) and primitive (M(ss), M(nn)) 
mesons. In this formula, and in other relations below, M's stand either for the mass 
or mass squared. 

The formula may be derived in the following way [H: For a meson nonet, the 
isosinglet qq mass matrix modified by the inclusion of the mixing amplitudes A, 
which are here assumed flavor-symmetric, for simplicity, 

^(2x2) _ ( M(ss) + A V2A \ 

M ~ { V2A M(nn) + 2A j ' 1 j 



1 In the pseudoscalar sector the axial anomaly and instantons can lift the degeneracy In the 
scalar sector only instantons lift the degeneracy, and in other sectors these effects are not significant 

1- 
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States 


I 


II 


III 


1- 


p, K*, uj, cj) 


1013.0 ±0.5 


1014.5 ±0.4 


1017.8 ±5.8 


2++ 


oa(1320), #|(1430), /a(1270), #(1525) 


1541.7 ±5.3 


1539.4 ±5.5 


1539.9 ± 10.6 


2-+ 


tt 2 (1670), ir 2 (1770), 772(1645), ?7 2 (1870) 


1870.7 ± 26.0 


1869.1 ±25.3 


1869 ±24 


3— 


#,(1690), # 3 *(1780), w 3 (1670), ^(1850) 


1863.1 ±9.2 


1862.3 ±9.0 


1863.2 ± 16.1 



Table I. Comparison of the values for M(ss) given by Eq. (1) for linear and squared 
masses (in columns I and II, respectively), and the quark model motivated relation 
discussed in the text for linear masses (in column III). All the mass values are given 
in MeV. All the mass values are taken from ref. ||. Uncertainties represent the 
combination of experimental uncertainties and the scale of electromagnetic splittings. 



is diagonalized by the masses of the physical mesons: Af(2x2) diag (M(ss), M(nn) 
Writing down the trace and determinant conditions for the matrix (2) and eliminating 
A from them leads to the formula (1). 

In Table I we compare the predictions of Schwinger's formula for the mass of the 
primitive ss state, based on the masses of the primitive nn state assumed to be the 
mass of the isovector, and the two physical states, with the predictions of the quark 
model motivated linear mass relation M(nn) + M(ss) = 2M(sn), where we use the 
mass of the isodoublet for M(sn). 

Schwinger's formula reflects the limitations on changes of masses from primitives 
to physical states that can occur only by mixing. To the extent that one may identify 
masses of primitives, by extracting mixing angles or using model assumptions (as 
above), the formula agrees well with data in most cases. 

We interject a comment on the primitives. They are not simply the result of 
suppressing quark loop contributions. For example, in quenched lattice QCD, OZI 
rule violating decays of ss states into two light quark- ant iquark pairs may still occur. 
Perhaps a useful analogy may be found in the issue of CP- violating kaon decays - 
there are both mass mixing induced and intrinsic CP-violating decays. Thus here 
it is insufficient simply to extract the mixing angle to the primitives by determining 
from the experimental decay rates which combinations of the physical states would 
not violate the OZI rule. One must allow for the analog of the intrinsic rule violating 
decays. 

Because of this, it continues to be much easier to define primitives on the basis 
of some model assumption, such as that the light isoscalar primitive is equal in mass 
to the light isovector state. Further, the notion that the light-strange mixed quark 
primitive averages the masses of the pure light and pure strange primitives receives 
adequate support from a wide range of model calculations dating back to the earliest 
days of the quark model through to the present refined potential and string models. 
The success of the Schwinger formula with these definitions of the primitives estab- 
lishes the existence of primitives with these properties. At this stage the glueball is 
"integrated out" in the formula. 

The agreement of Schwinger's formula with data is, however, problematical, as 
the mixing must pass through glueballs determined by QCD dynamics. Here we shall 
examine such mass relations by explicitly including the nearest glueball intermediate 
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states to describe the mixing and complete the phenomenological picture, noting that 
when the glueballs are near to the mesons one should not be able to "integrate out" 
their effects. We show that if the single intermediate glueball state is only modestly 
more massive than mesons in the light sector, then the 2x2 and 3x3 matrices 
give almost identical results. This is to be expected in the large glueball mass limit, 
but we find it develops surprisingly early - just above 2 GeV. We describe how the 
situation generalizes in the presence of multiple meson nonets. 

The validity of Schwinger's original formula (1) should be expected to break down 
in cases when additional isosinglets mix strongly with the primitive mesons. Such 
additional isosinglets may be glueballs or hidden flavor heavy quark mesons with the 
quantum numbers of the nonet mesons. In such cases, one has to consider a larger 
mass matrix, N x N, N > 2, for which an analog of the original formula (1) may be 
expected to exist. 

In this paper we derive such an analog of Schwinger's original mass formula (1) 
in the case of any number of additional isosinglets (we need not specify their nature 
for our present purposes). 



2 Generalized Schwinger mass formula 

Let us first note that a field theory in general, and QCD in particular, does not 
have a mass matrix corresponding to composite states. Indeed, the diagonalization of 
a mass matrix corresponds to a linear transformation of fields that diagonalizes the 
Lagrangian, which can only be done for Lagrangians which are quadratic in fields, e.g., 
the Klein-Gordon Lagrangian, with quadratic interactions. A mass matrix appears 
therefore only in effective field theories of this type. However, precisely this type of an 
effective meson theory may be expected to exist for QCD. In fact, the emergence of 
such a theory was shown in |7|] through the bosonization of QCD. Hence, we assume 
the existence of an effective theory that describes the mixing of nonet mesons with 
additional isosinglets. 

We start with the case of one additional isosinglet. Assume, for simplicity, that 
an additional isosinglet (glueball or hidden flavor heavy quark meson), i, couples to a 
pair of primitive (nonet) mesons. Then, we have the following 3x3 hermitian mass 
matrix where, in addition to the quark mixing amplitudes, we have the amplitudes 
of the i — qq mixing (the couplings) which we denote by B : 



^(3x3) = B * M ^ +A Ar (3) 




where M(i) is the mass of the primitive i. From now on, in contrast to Eq. (2), we 
allow for possible violation of exact flavor SU(3) symmetry in terms of arbitrary r. 
The SU(3) limit will then correspond to r = The form of this mass matrix is 
motivated in ref. ||. The mass matrix (3) was also discussed in |§, and its SU(3) 
symmetric version in fT0|j . 
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Note that in Eq. (3), as well as in Eq. (2), A must be real. With B and r both 
real the matrix (3) is the most general parametrization of 3 x 3 (real) symmetric 
matrix, since it contains 6 independent parameters. 

The mass matrix (3) is diagonalized by the masses of the three physical states, 
M {3x3) diag (M(i), M(ss), M(nn)) , which are determined from the three eigen- 
value equations (which follow from Det (^M.^ 3x3 "> — A/) = 0). Eliminating A and B 
from the eigenvalue equations leads, upon some algebra, to the formula 



1 + \r\ 2 )M(i) - \r\ 2 M{ss) - M(nn) (1 + \r\ 2 )M(ss) - \r\ 2 M(ss) - M(nn) 



(1 + \r\ 2 )M(nn) - \r\ 2 M(ss) - M{nn) 



(1 + \r\ 2 )M(i) - \r\ 2 M(ss) - M(nn) 



M(ss) - M(nn) 



0. 



(4) 



which is the generalized Schwinger mass formula for the case of a meson nonet with 
an additional isosinglet. We note that the new formula does not depend on either the 
values of A or B, just as the original formula does not depend on the value of A. 
Under the conditions that either 

(i) (1 + \r\ 2 )M(i) > \r\ 2 M(ss) + M(nn) (see Eq. (4)), or 

(ii) \B\ < A, 

so that M(i) ~ M(i), the generalized Schwinger formula (4) (with r = y2) reduces 
to the original one, Eq. (1). 

In the case (iii) A = 0, 
i.e., no direct coupling between mesons, one can combine the new formula (4) with 
one of the eigenvalue equations (the trace condition for the matrix (3)) 



M(i) + M(ss) + M{nn) = M(i) + M(ss) + M{nn) 



(5) 



in order to determine two unknown masses. Under the assumption that (iii) is the 
realistic case for meson spectroscopy, i.e., that the observed meson mass spectra result 
from coupling of the primitive qq states through an additional isosinglet, one can 
calculate the masses of two physical states, e.g., M(i), M(ss), as functions of M(i), 
using the known M(nn), M(ss). In Figs. 1,2 we compare M(ss) calculated in this 
way with the range allowed for M(ss) by experimental data 0, for the J PC = 3 
multiple^ (h ere we use Eqs. (4), (5) for linear masses). One sees that the value 
calculated lies within the range allowed by data for a rather large bare mass of an 
additional isosing let, ~ 2 GeV. 

Note that all isosinglet masses above 2.5 and 2 GeV, respectively, are consistent 
with M(ss) from experiment, so that the mass of an additional isosinglet cannot be 
uniquely determined from current experimental data. Thus, the procedure described 
above (only) fixes the lower bound on the mass of an additional state (i). 

2 The other meson nonets can be analyzed in a similar way. The case of the J PC — ++ has been 
studied in detail in ref. |Ol. 
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3 Further generalization of the new Schwinger mass 
formula 



The most general case for a mass formula of the type (4) is evidently the case of the 
mixing of arbitrary numbers of both mesons and additional isosinglets with common 
J PC . This most general case can be decomposed into the following two cases: (i) an 
arbitrary number of additional isosinglets (glueballs and hidden flavor heavy quark 
mesons) mixing with one pair of (nonet) mesons built from u, d, s quarks, and (ii) an 
arbitrary number of pairs of excited mesons (with common J PC ) built from u, d, s 
quarks mixing with one isosinglet (glueball or hidden flavor heavy quark meson) . We 
now examine these two cases in more detail. 



3.1 The case of an arbitrary number of additional isosinglets 



In this case the new Schwinger mass formula has an explicit further generalization, 
as follows. 

Consider the mixing of (N—2) extra isosinglets with a pair of nonet mesons and in- 
troduce, as previously, the corresponding bare masses M(ii), . . . , M(i]y- 2 ), couplings 
Bi, . . . ,£>tv-2, quark mixing amplitudes A and flavor SU(3) breaking parameter r. 
The N x N, N > 3 hermitian mass matrix is then diagonalized by the masses of the 
N physical states: 



/ M(ii) 






B* 
\ B{r* 





M(i 2 ) 





B* 
B* 2 r* 








M(i N 
B 



Br 
B 2 



B ir 
B 2 r 



\ 



B 



N-2 
N-2 r 



B N _ 2 B N _ 2 r 
M(ss) + A Ar 

Ar* M{nn) + A\rf 



(6) 



diag (M(ii), M(i 2 ), M(i N _ 2 ), M(ss), M(nn)) . 



(7) 



Here we assumed that it is possible to choose a basis where the primitive additional 
isosinglet states are already diagonalized. For example, if the three isosinglets are 
a glueball, cc and bb, such a basis can to a good approximation be chosen due to 
the small mixing between the glueball and hidden flavour heavy quark mesons. We 
also assumed that the nn meson to ss meson coupling ratio of all isosinglets is a 
universal constant r. In the limit of SU(3) symmetry this assumption becomes exact, 
and r = \/2. 

Since the initial N x N mass matrix (6) can be rewritten in the following effective 
3x3 form: 

M(i) B Br 

Bt M(ss) + A Ar | , (8) 

BV* Ar* M(nn) + A\rf 
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where 



M(i) = 



( m{h) o 

M(i 2 ) 



V o 








M(ijv_ 2 ) / 



B 



B 2 



(9) 



it is evident that the general NxN case effectively reduces to the 3x3 case considered 
above. Applying the methods of linear algebra, one can obtain upon some calculation 
the following most general Schwinger type mass relation for the NxN case:0 



Af-2 



I [ ((1 + \r\ 2 )M(i k ) - \r\ 2 M(ss) - M(nn))({l + \r\ 2 )M(ss) - \r\ 2 M(ss) - M(nn 

k=l 

x ((1 + \r\ 2 )M(nn) - \r\ 2 M(ss) - M(nn)) 

Af-2 2 

+ [] U 1 + M 2 )M{i k ) - \r\ 2 M{ss) - M(nn)) \r\ 2 (m(ss) - M(nn)) =0. (10) 

k=l 

In the case of only one additional isosinglet, Eq. (9) reduces to the generalized 
Schwinger mass formula (4). 



4 The case of an arbitrary number of mesons 

As discussed above, the generalized Schwinger mass formula (10) reduces to the orig- 
inal formula (1) in the limits where the primitive mass of an additional isosinglet 
goes to infinity, or the isosinglet-meson coupling goes to zero. The main use of the 
new Schwinger formula is hence where an additional isosinglet is in the vicinity of 
the mesons, where the conventional Schwinger formula is expected to be inapplicable. 
However, as the two examples of the J PC = 2 ++ and 3 multiplets considered above 
show, the mass of such an additional isosinglet must generally be expected to be 



large (in agreement with, e.g., lattice QCD calculations of glueballs |TT[). Hence the 
mesons in the new Schwinger formula would need to be radially or orbitally excited 
mesons, or hybrid mesons (in all nonets except the scalar nonet (T^|). It now becomes 
interesting to consider whether the excited states could influence the mixing of an 
additional isosinglet with the ground states. It should be decided how many pairs 
of excited states to incorporate. There are arguments that the effective number of 
hadronic states of any given type is limited El . We shall discuss a 7-state case which 



is sufficiently complicated to simulate the general case. The 7x7 matrix has also 
been discussed in a different context in ref. 



Here we investigate a set of three pairs of mesons and test for the effect by choosing 
an additional isosinglet at various masses. We choose ground states at M(nn) = 1.2 
GeV, M(ss) = 1.5 GeV, first excited states at M«(nn) = 1.7 GeV, M^(ss) = 2 

3 One of us (L.B.) wishes to thank CD. Levermore for a very useful discussion on this point. 
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GeV, and second excited states at M {2 \nn) = 2.1 GeV, M^(ss) = 2.4 GeV. The 
mass splittings chosen are typical of constituent quark models 0. The additional 
isosinglet mass cases we consider are (i)-(vi). Each of the following two cases, 

(i) M{i) = 1.35 GeV, 

(ii) M{i) = 1.6 GeV 

simulates the scalar meson case for which an additional isosinglet (the scalar glueball) 
is in the vicinity of the ground state mesons [12 1 . Each of the following three cases, 
(hi) M(z) = 1.85 GeV, 

(iv) M(z) = 2.05 GeV, 

(v) M(i) = 2.25 GeV 

simulates the tensor meson case for which the mass splitting between the tensor 
glueball and the nn ground state meson is [TT| ~ 0.9 GeV. Finally, the case 

(vi) M{i) = 2.5 GeV 

simulates the 2 h and 3 cases for which the mass splitting between the correspond- 



ing glueball and nn ground state meson is W% ~ 1.3 GeV. 



For each of the above six choices of the primitive mass of an additional isosinglet, 
we first consider the full 7x7 mass matrix for the mixing of the isosinglet with 
ground states and both pairs of excited states which generalizes the 3x3 matrix 
(5). The couplings of an additional isosinglet to mesons are taken to be the same 
for ground states and both excited states, for simplicity:^ B = 0.4 GeV 2 , and no 
direct meson-meson coupling^: A = 0. We then consider separately the mixing of an 
additional isosinglet with ground states and first and second excited states, each of 
which is described by the corresponding 3x3 mass matrix with the same B = 0.4 
GeV 2 and A = 0. In each case, we calculate the masses of all the physical states, 
using the appropriate quadratic mass relations. Our results in the full 7x7 case and 
the separate 3x3 cases are compared in Table II. 

As seen in Table II, the largest discrepancy between the results of the 7x7 and 
3x3 cases occurs for the nn ground state in the case (i) and constitutes ~ 8%. 
In every other case, the corresponding discrepancy does not exceed ~ 2% (which 
corresponds to a mass uncertainty of < 50 MeV), which means that the generalized 
Schwinger mass formula in every 3x3 sub-basis should be very accurate. We note 
that in the actual case of the scalar nonet with smaller B = 0.3 GeV 2 HT2T1, in the case 



(i) one would obtain M(nn) = 1032 MeV, M(ss) = 1535 MeV for the 7 x 7 case, and 
M(nn) = 1068 MeV, M(ss) = 1565 MeV for the 3x3 case, so that the discrepancy 
between the results for the M{nn) reduces to only ~ 3.5%. 

We have also done a similar comparison between the 7x7 and 3x3 cases for 
the same six choices of the primitive mass of an additional isosinglet and the same 
bare mesons masses, but with nonzero direct meson-meson couplings, A, in each of 
the 2x2 mesons sub-bases, in addition to nonzero B. For A — B — 0.4 GeV 2 , the 
largest discrepancy between the corresponding results again occurs for the nn ground 
state in the case (i) and constitutes ~ 7%. In every other case, the corresponding 
discrepancy does not exceed ~ 5% (which corresponds to a mass uncertainty of < 80 

4 This value represents the realistic value of the glucball-mcson coupling, as extracted from data 



by the phenomcnological analysis in ref. 15 . 

5 It can be shown that there are analytical Schwinger-type formulae in this case pi 
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MeV). 

Note that in the above analyses, for ground states and both pairs of excited states 
we formally used a common value of B. However, the couplings of the excited states 
(to the glueball) may be smaller. This would give rise to even smaller errors than 
quoted above. 

On the basis of both possibilities for zero and nonzero A analyzed, we conclude 
that the generalized Schwinger formula is only weakly affected by the inclusion of the 
excited states, and must therefore hold with good accuracy for almost every 3x3 
isosinglet-mesons mass sub-matrix. 



5 Conclusions 

Generalized Schwinger mass formulas have been derived for an arbitrary number of 
isosinglets coupled to a meson nonet. The formulae have significant utility as they 
mainly refer to experimentally measurable masses, and not unknown couplings; they 
can therefore be used to predict unknown masses. For real couplings, the generalized 
Schwinger mass formula with the inclusion of a single isosinglet describes the most 
general parametrization of a symmetric 3x3 matrix. In other cases, the matrices are 
not the most general hermitian matrices, albeit physically relevant ones. 

When multiple meson excitations were included in addition to a single isosinglet, 
we considered the most general matrices. It has been shown for some numerical 
examples that to a good approximation it is sufficient to restrict to the relevant 3x3 
isosinglet-mesons subspace. 
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Figure 1 
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Fig. 1. Comparison of the mass of the physical ss state calculated from Eqs. (4), (5) 
as a function of the bare mass of an additional isosinglet, M(i), with the range allowed 
by experimental data of ref. (horizontal band), for the J PC = 3 multiplet. 
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Figure 2 
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Fig. 2. The same as Fig. 1 for the J PC = 3~ multiplet. 
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case (i) 


case (ii) 


case (iii) 


case (iv) 


case (v) 


case (vi) 




7x7 


3x3 


7x7 


3x3 


7x7 


3x3 


7x7 


3x3 


7x7 


3x3 


7x7 


3x3 


M(nn) 


918 


994 


1051 


1086 


1115 


1131 


1141 


1150 


1157 


1162 


1169 


1171 


M(ss) 


1545 


1607 


1405 


1444 


1452 


1468 


1470 


1478 


1479 


1484 


1486 


1488 


M^\nn) 


1765 


1763 


1802 


1802 


1617 


1581 


1645 


1632 


1662 


1657 


1674 


1671 


M^(ss) 


2015 


2020 


2019 


2030 


2028 


2065 


1929 


1947 


1964 


1973 


1981 


1985 


M^\nn) 


2134 


2127 


2148 


2135 


2179 


2155 


2235 


2195 


2058 


2015 


2073 


2060 




2410 


2409 


2412 


2411 


2418 


2415 


2428 


2423 


2463 


2446 


2373 


2364 


M(i) 


1344 


1342 ±51 


1574 


1631 ± 96 


1870 


1841 ± 76 


2042 


2007 ±87 


2317 


2292 ± 14 


2603 


2544 ± 23 



